Intercalation and buckling instability of DNA linker within locked chromatin fiber 
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The chromatin fiber is a complex of DNA and specific proteins called histones forming the 
first structural level of organization of eukaryotic chromosomes. In tightly organized chromatin 
fibers, the short segments of naked DNA linking the nucleosomes are strongly end constrained. 
Longitudinal thermal fluctuations in these linkers allow intercalative mode of protein binding. We 
show that mechanical constraints generated in the first stage of the binding process induce linker 
DNA buckling; buckling in turn modifies the binding energies and activation barriers and creates a 
force of decondensation at the chromatin fiber level. The unique structure and properties of DNA 
thus yield a novel physical mechanism of buckling instability that might play a key role in the 
regulation of gene expression. 
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Notwithstanding its evident biological importance, 
DNA is a fascinating object for physicists due to its re- 
markable physical properties. First, DNA is a molec- 
ular spring exhibiting stretch, twist, and bend elastici- 
ties. These elastic properties have been thoroughly stud- 
ied both theoretically [0 and by single molecule experi- 
ments that led to unexpected results as for instance the 
structural transitions observed when a sufficiently strong 
pulling force is applied to the DNA molecule Sec- 
ondly, DNA is a polyelectrolyte and is involved in nu- 
merous electrostatic effects ||; in particular, its nega- 
tive charge density is enough large to attract a sheath 
of counterions almost independent of the salt concentra- 
tion, a phenomenon known as Manning's condensation 
j|. Third, the DNA double helix undergoes a denat- 
uration transition studied both theoretically || |q] and 
experimentally 0. We here focus on a fourth specific 
feature of DNA, namely, intercalation ||, that allows the 
binding of planar molecules between adjacent DNA base 
pairs (see Fig. 1). 




Figure 1: fntercalative mode of binding within DNA for 
a mono-intercalator (left) and a bis-intercalator (right). 
One (respectively two) domain(s) of the binding protein 
comes in between two successive base pairs, inducing a 
rise and unwinding — At. For instance, ethidium bro- 
mide is a mono-intercalator with A^ = 2A and At = 26°. 



Quite recently, theoretical Q [jl(| (jTl| and experimen- 
tal groups |T^] |13| started to explore the nanomechanics 
of chromatin since free DNA is not the relevant instance 
of DNA in vivo. In the nuclei of plant and animal cells, 
DNA is actually organized in a hierarchy of structural 
levels. The first one is the wrapping of 146 bp of DNA 
around histone octamers to form the nucleosomes. Nu- 
cleosomes remain connected by naked DNA segments — 
the linkers — of length between 10 bp and 100 bp accord- 
ing to the species and cell type. The nucleosome-dressed 
DNA molecule is further organized into a helical folding 
of about 30 nm in diameter that is called "the chromatin 
fiber" H] . We recently proposed on mechanical grounds 
that the chromatin fiber might exhibit a columnar 
packing of nucleosomes similar to columns observed in 



colloidal solutions of mononucleosomes |15(] . It suggests 
that chromatin might be locked into a strongly organized 
structure, induced by interactions within stacked nucle- 
osomes and secured by histone tails (Fig. 2). We call 
"locked chromatin" such a structure in which the ends 
of each linker are fixed in space due to strong three- 
dimensional positioning of nucleosomes within the fiber. 
A quite similar structure has been suggested long ago by 
Worcel et al. |l6| . We claim that it provides a plausible 
structure for facultative heterochromatinc. 




Figure 2: Proposed model of condensed chromatin, 
locked by interactions between stacked nucleosomes and 
presumably also by histone tails [^) . The chromatin fiber 
axis is here vertical. 



In this paper, we propose and describe a different 
mechanical property occuring specifically within con- 
strained DNA as encountered in locked chromatin. It 
consists of a buckling instability generated by longitudi- 
nal thermal fluctuations stabilized by intercalation. 

Intercalative mode of binding plays an important role 
in vivo; for instance the TBP (TATA-box-binding pro- 
tein) , which binds on specific sequences called the TATA 
boxes and plays a seemingly universal role in eukaryotic 
transcription initiation, is a (multiple) intercalator Jl7| . 
Intercalation is additionally involved in experimental in 
vitro or in vivo studies through the use of fluorescent 
dyes, as for instance ethidium bromide In the classi- 
cal model of intercalation, the binding process is decom- 
posed into three steps [Q p9| . The first step is a ther- 
mally activated local opening of DNA, creating a binding 
site for an intercalating molecule. The activation barrier 
is of order AG° con ^ w 6.5 kT according to Chaires |l^] 
|l9| . This opening is achieved through a local stretching 
Al > of the interbase-pair distance and a local un- 
winding —At (Fig. 3). More precisely, it corresponds to 
a two-state local conformational transition of the binding 
site involving changes in the DNA backbone and orien- 
tation of the bases and sugars pOf. The second step is 
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the insertion of intercalator inside the binding site (hy- 
drophobic transfer process) followed in a third step by 
the formation of molecular interactions, namely hydro- 
gen bonds between the intercalator and surrounding base 
pairs. The intercalated site is then slightly smaller in 
size than the preintercalation site required to first acco- 
modate the intercalator. We underline that the binding 
site opening is not induced by the presence of intercala- 
tor but merely by thermal fluctuations, then stabilized by 
intercalator insertion: the mechanism is not an induced 
fit but rather a conformational capture, in which interca- 
lation captures an "excited state" reached spontaneously 
by thermal fluctuations pi). 



5 kT 




Figure 3: Schematic drawing showing the modification 
of intercalation energy barrier and energy well coming 
from mechanical contribution AAG. The units are real- 
istic: the barrier is of order 7 kT; the excess energy is 
there about 2 — 3 kT. The binding energy is of order 
—20 kT; the excess energy is there about 1 — 2 kT (the 
intercalated site is smaller than the "open" site required 
in the first step of intercalation, which reduces AAG into 
AAG mt ) @ M. 



Let us now consider an intercalation within a linker 
embedded in a locked chromatin fiber, hence with trans- 
lationally and rotationally fixed ends. We describe the 
linker DNA within a standard continuous model (gen- 
eralized worm like chain) as an extensible rod of length 
I whose mechanical properties are fully described by its 
bending persistence length A ss 50 nm, its twist persis- 
tence length G « 75 nm and its stretch modulus 7 w 1200 



pN jy. Since the distance between linker ends is fixed, 
prescribed by the global architecture of the chromatin 
fiber, the first step of intercalation generates mechanical 
constraints in the remaining part of the linker, namely, a 
compression — Al and an overtwist Ar, from which fol- 
lows an excess energy AAG. The formation of an inter- 
calation site thus requires to overcome the energy barrier 
AG con/ = AG° conf + AAG, with 



AAG = 7 



(AQ 2 
2(1 - h) 



kTC 



(Ar) 2 
2(1 - h) 



(1) 



where h is the interbase-pair distance. We estimate this 
excess energy with h = 3.4 A , A^ = 3 A , At = 30° (we 
take these structural values slightly larger than those ob- 
served after intercalation, presuming that the preinter- 
calation site should be larger in size). We find AAG « 
2.4 kT for I = 30 bp (10 nm) and AAG w 3.5 kT for I = 
20 bp (6.7 nm). This excess energy is of order of thermal 
energy provided I is not too small; it shows that ther- 
mally activated creation of an intercalation site within 
a constrained linker is still possible for I greater than a 
minimum length « 20 bp. We underline that the ther- 
mal fluctuations here invoked are longitudinal fluctua- 
tions; bend fluctuations play a negligible role due to the 
short length of the linker compared to the bend persis- 
tence length of DNA |2^] . Mechanical constraints do not 
significantly modify the energy terms involved in the fol- 
lowing steps, namely, the free energy cost associated with 
the hydrophobic transfer step, the energy associated with 
modification of counterion surrounding and the energy of 
newly established molecular interactions between interca- 
lator and DNA (l8) . Note that after intercalation, AAG 
is to be computed with final values of Al and Ar, which 
gives a "mechanical" correction AAGint to AG that is 
smaller than AAG. Mechanical constraints modify both 
the thermodynamics of the binding (depth of the energy 
well) and the kinetics of the binding (height of the ac- 
tivation barrier): the activation barrier is increased by 
an amount AAG whereas the binding energy rises by a 
smaller amount AAGi„ t (see Fig. 3). 

Compared to intercalation within free DNA, the differ- 
ence lies also in the stresses generated when an interca- 
lation site opens in an end-constrained linker. The force 
induced by the compression — Al of the linker outside the 
intercalation site writes 



F int = 7 



Al 



(2) 



For I = 30 bp. 
F mt = 36 pN. 



i.e. / = 10 nm, and Al — 3 A , we obtain 



It is well known that above a critical threshold, the 
compression of a column shifts into bending, what is 
called "buckling." This phenomenon persists up to 
nanometer scale p3| . In order to determine whether the 
force Fi n t will induce, or not, buckling of the linker, we 
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compute the critical force F c corresponding to buckling 
instability threshold. At this stage, the boundary con- 
ditions at the linker ends, namely, the mode of anchor- 
ing of linkers onto the nucleosomes, have to be described 
pre cisely. Standard computation of material mechanics 
p4f shows that boundary conditions might be taken into 
account through a numerical coefficient v reflecting the 
degree of end fixing, e.g. v = 1/2 for clamped ends and 
v = 1 for hinged ends; any value of v > 1/2 is in fact 
possible for various elastic joints. The critical Euler load 
writes accordingly 

TT 2 kTA 



F„ 



H) 2 



(3) 



where A is the above-mentionncd DNA-bending persis- 
tence length. Buckling takes place if Fi nt > F c or equiv- 
alently if 



I > 



K 2 kTA 
v 2 jAl 



h 



ir 2 kTA 
v 2 1 Al 



= I, 



(4) 



This criterion could have been obtained by comparing 
the compression energy j(Al) 2 /2(l — h) with the critical 
buckling energy lF 2 /2^. We get 



for v = 1 F c w 20 pN and l c « 16 bp 
for v = 1/2 F c « 80 pN and l c « 65 bp 



(5) 



hence intercalation induces buckling in typical linkers (30 
bp) if their ends are hinged but not if they are clamped. 
This result shows that for linker length values between 
30 and 60 bp, buckling may be selectively controlled by 
the anchoring of the linkers at the entry-exit points on 
the nucleosomes. Indeed, according to the biochemical 
status of the nucleosome core, either linker DNA is tightly 
grafted onto the core, which corresponds to clamped ends 
(y = 1/2), either DNA might unwrap, typically by 5- 
10 bp on each side, at negligible energetic cost 13 [p5| , 
which corresponds to hinged ends (y = 1) (Fig. 4). 





Figure 4: Buckling of the linker is possible only when 
the linker ends are allowed to unwrap from the nucleo- 
somes. Unwrapping by 5 bp on each side is enough to 
accomodate a buckling angle a — 20°, associated with 
Al/l = 0.03; within these bounds, the situation is me- 
chanically equivalent to hinged ends {y = 1). 



Linker buckling in turn modifies the energetics of in- 
tercalation. When buckling occurs, excess energy is 



now the sum of three contributions: the compression 
energy before buckling j(Alb) 2 j2(l — h), where Alt = 
F c (l — h)/"f is the stretching deformation required to 
reach the buckling threshold, the compression energy af- 
ter buckling F c (Al — Alb) and the overall twist energy 
kTC(Ar) 2 /2(l - h). It comes 



AAG 



buckled 



= F n 



Al 



FS - h) 
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kTC(Ar) 2 
2(1 - h) 



(6) 



This expression points out that once a linker is buck- 
led, the energy cost to stretch it further is linear in 
the length increase and proportional to 1/v 2 , hence 
tunable by a modification of the end status. Buck- 
ling thus facilitates the insertion of additional interca- 
lators, and perhaps more crucially, it allows the inser- 
tion of multi-intercalators as, for instance, the above- 
mcntionned TBP. 

At the fiber level, the buckling force F c induces 
stresses, that are easy to estimate within our structural 
modeling of the condensed chromatin fiber. When the 
two linker ends are anchored in the same way on the nu- 
cleosome core, the force acts along the intercalated linker; 
it exhibits a component F c cos z along the chromatin fiber 
axis, where z is the angle between the linkers and this axis 
H . Buckling also generates a radial shear whose exact 
expression depends on the details of the chromatin fiber 
structure. In the proposed locked structure of the fiber 
(Fig. 2), cosz is around 0.6, close to its maximum value: 
this structure is thus at the same time the more easy to 
lock thanks to interactions between stacked nucleosomes 
and the more easy to open thanks to the decondensing 
force generated by intercalation-induced linker buckling. 
This two-fold property provides additional support of the 
biological relevance of the proposed structure. Note that 
once the linker is buckled, the force experienced by the 
linker, hence the force of decondensation does not vary 
significantly with the number of bound intercalator pro- 
teins. 

Buckling and decondensing forces are generated by the 
hyperstatic structure of the locked chromatin fiber. Ac- 
tually, intercalation induces not only compression but 
also torsional strains within the linkers (which hence be- 
have not only as columns but also as shafts). Twist- 
induced stresses are yet proportional to the number of 
intercalated molecules and can produce a decondensing 
force of several piconewtons along the chromatin fiber 
axis, but buckling, occuring before multiple intercalation, 
modifies the picture and bypasses twist effects by turning 
them into writhe. We checked that the resulting contri- 
bution plays a secondary role. 

Taken together, the above results suggest the following 
scenario providing the first physical bases for the decon- 
densation process: 

(i) Thermal fluctuations allow the formation of an in- 
tercalation cavity which induces a compression force in 
the remaining part of the linker. According to the mode 
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of linker anchoring onto the nucleosome core, the com- 
pression force may induce buckling of the linker. 

(ii) Buckling is further stabilized by insertion of an 
intercalator into the cavity. Moreover, buckling allows 
the insertion of further intercalators in the same linker 
so that the limiting step is indeed buckling. 

(iii) The buckling force generates in turn decondensa- 
tion forces at the chromatin fiber level: a stretching force 
along the axis of the chromatin fiber as well as shearing 
forces between stacked nucleosomes. 

This nonlinear decondensation mechanism is irreversible; 
chromatin condensation will occur along another path- 
way, involving electrostatics to induce compaction of the 
fiber (currently under study). 

The physical mechanism here proposed might be of 
biological importance in the regulation of transcription. 
Transcription initiation requires a step of decondensa- 
tion of the chromatin fiber in order to give access to 
transcription machinery p6| . In this paper, we raised 
the point that binding of transcription factors and other 
DNA-binding proteins that monitor the early stages of 
transcription should occur within a locked fiber, hence 
in mechanically constrained linkers. Structural modifica- 
tions of linker DNA here play a role through the forces 
they generate, without any ATP-consuming mechanism. 
This point has yet been put forward to account for the 
cooperativity of protein binding on naked DNA j27J. 

Presumably, the chromatin fiber architecture has been 
deviced in the course of evolution not only to efficiently 
pack DNA inside the nucleus, but also, as proposed here, 
to play a mechanically active role in gene expression reg- 
ulation. 
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